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3. 


ABSTIACT 

Progress  made  during  ths  second  quarter  of  a  research 
program  directed  toward  a  fundamental  study  of  the  non-steady 
ccabustion  of  solid  propellants  with  aiq>lioation  to  rocket  instab¬ 
ility  is  described. 

A  film  strip  of  a  metallised  solid  propellant  burning 
in  an  oscillatdbog  pressure  field  was  analyzed-in  order  to  determine 
the  phase  relationship  between  the  inposed  pressure  maxima  and  an 
apparent  %fave  of  luminosity  in  the  co^ustion  gases.  The  inplications 
of  these  observations  as  they  affect  the  model  of  the  burning  propel¬ 
lant  are  briefly  considered. 

Two  types  of  experimental  observation  which  will  be  used 
to  study  non-steady  combustion  are  outlined.  One  method  involves 
the  measurement  of  the  temperature  of  the  cosbustion  gases  down¬ 
stream  from  the  flame  zone.  The  use  of  a  "particle  track"  method 
for  the  observation  of  pressure-velocity  relationships  in  the  burnt 
gases  is  also  considered. 
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1.  IKTRODUCTION 

Tho  complex  interaction  of  on  oocilloting  preosuro  field 
with  tho  flame  of  a  burning  solid  propollont«  which  opparcntly 
lies  at  tho  heart  of  tho  non-stcody  burning  problem,  is  given 
further  consideration  in  this  second  quarterly  progress  report. 
Earlier  work,  in  which  tho  feasibility  ond  efficacy  of  a  number  of 
frequently  discussed  experimental  techniques  wore  investigated, 
highlighted  tlio  difficulty  of  contriving  oxporimonts  capoblo  of 
yielding  sufficiently  sensitive  mcasurcri^nts  of  tho  acoustic  re¬ 
sponse  of  a  burning  surface.  In  this  report,  experiments  are  de¬ 
scribed  which  seem  capable  of  yielding  results  pertinent  to  the 
acoustic  response  problem.  At  tho  same  time,  experimental  data 
obtained  in  a  separate  investigation  have  boon  reduced  in  a  new 
way  to  give  information  regarding  tho  phase  relationship  between 
a  pressure  oscillation  and  a  response  of  tho  cembustion  process. 


ZI.  FILM  Al^'ALYSZS 

The  burning  of  metals  incorporated  in  a  solid  propellant 
matrix  has  been  under  investigation  by  Dr.  William  A.  Wood  of  Rohm 
&  Haas  Co.  (Redstone  Division)  for  some  time.  As  part  of  this  study, 
cinematography  v/as  employed  to  record  tho  burning  process  of  a 
metal-containing  propellant  under  conditions  of  oscillating  chainber 
pressure.  Tho  experimental  setup  is  as  shown  in  Figure  1.  It 
^s  noted  that,  in  tho  combustion  gases  from  a  propellant  containing 
TA  magnesium,  there  were  clearly  visible  luminous  zones.  Figure  2 
«  a  sketch  of  a  typical  framo^  showing  the  dark  and  luminous  zones. 
These  zones  moved  with  the  product  gases,  but  tho  separation 
between  them  oscillated  with  time  at  the  same  frequency  as  the 
pressure  oscillations.  Tho  experimental  conditions  for  the  par¬ 
ticular  film  strip  studied  are  outlined  in  Table  1. 

The  phase  difference  between  the  luminous  waves  and  the 
pressure  cycle  is  of  fundamental  importance  to  an  understanding 
of  the  combustion  mechanism.  Unfortunately,  no  pressure-time 
history  was  available  with  the  film  strip.  However,  it  was  pos¬ 
sible  to  obtain  the  pressure-time  history  in  the  cavity  from  the 
variation  in  the  spacing  of  the  luminous  zones.  These  variations 
in  spacing  arise  from  oscillations  in  the  gas  density  which  are 
in  turn  related  to  the  pressure  oscillations  in  the  gas  column. 

An  analysis  was  made  of  ^d}o^t  500  frames,  from  which  was  found  an 
average  phase  lag  of  130°  between  a  pressure  maximum  and  an  emis¬ 
sion  of  a  luminous  zone  from  the  burning  surface. 

In  general  a  phase  lag  might  be  caused  by  at  least  two 
conditions.  First,  an  oscillatory  variation  in  the  propellant 
combustion  process  might  result  in  the  emission  of  a  type  of  ''smoke" 
periodically.  Of  course,  periodic  increases  in  pressure  cause 
corresponding  increases  in  gas  evolution  rate  at  the  surface, 
but  it  is  the  variations  in  mass  fraction  of  smoke  that  would 
be  carried  with  the  stream  and  show  up  in  the  photograph.  In 
this  case,  the  MgO  particles  could  become  visible  by  acting  as 
reflectors  of  the  light  fran  the  burning  propellant  surface.  The 
phase  relation  here  can  only  be  conjectured. 
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Second*  the  luminous  zones  could  represent  temperature 
peaks  resulting  from  entropy  v/aves.  These  entropy  waves  are 
carried  with  the  stream  and  are  caused  by  tl»e  fact  that  the  flame 
temperature  holds  constant  while  the  pressure  oscillates.  If  so. 
the  phase  lag  is  expected  to  be  180O.  The  finding  of  a  phase 
lag  of  130®  is  obviously  a  matter  to  be  investigated  in  farther 
e3q)erlraents.  These  experiments  would  distinguish  between  thermal 
luminosity  or  "smoke"  as  the  fundamental  cause  for  the  visible 


ZZZ.  THE  ACOUSTZe  ADMITTANCE 


McClure  and  coworkers  have  suggested  that  the 
acoustic  admittance  "Y"  of  the  burning  surface  is  an  important 
parameter  in  determlng  the  stability  of  a  rocket  chamber.  The 
equation  for  acoustic  admittance  may  be  derived  as  follows  t 

»  /O/y  (II 


.  A  d  nr 
^  P  dp 


+ 


Where t  m  =  mass  flow  rate  per  unit  area 
/O  “  gas  density 
V  •  velocity 
p  B  pressure 

But  "Y",  the  acoustic  admittance  is  defined  byi 


hence; 


If  the  period  of  the  oscillation  is  long  compared  to  the  combustion 
time  of  the  reactants*  the  combustion  gas  will  be  emitted  from  the 
flame  zone  with  a  temperature  which  is  the  same  at  all  times*  since 
the  adiabatic  flame  temperature  is  essentially  Independent  of  pres- 


For  a  constant 


flame  temperature* 

- 


6 


leading  to 


Jt^/m 

dp/p 


writing  yn.  a  and  d  • 

I® 

arrive  at  y*  -id  *  * 

P 


%m  imedlately 

(5; 


Equation  5  may  be  compared  withi 


ae  derived  by  McClure*. 

The  difference  is  that  McClure  took  the  Lagrangiaa 
preaeure-deneity  derivative  while  in  this  report  the  Bulerian 
derivative  was  used.  It  is  believed  that  this  approach  is  the 
correct  one. 

This  comparison  shows  that  the  use  of  a  constant  flame 
temperature  implies  greater  stability«  although  the  presence  of 
other  gain  and  loss  mechanisms  means  that  stability  or  instability 
of  a  rocket  chamber  ceui't  be  ea^lained  in  terms  of  the  acoustic 
admittance  of  the  propellant  alone. 


IV.  ENTROPY  mVE  HYPOTHESIS  AMD  SUGGESTED  BXFBRIMEMTS 

In  the  derivation  of  equation  (5)  it  was  assumed  that  ‘ 
the  combustion  gas  is  emitted  from  the  flame  at  essentially  the 
same  temperature  throughout  the  pressure  cycle.  However,  when 
the  gas  enters  the  duct,  its  temperature  will  oscillate  in  phase 
with  pressure  as  illustrated  in  Figure  3.  In  this  time-temperature 
plot  it  can  be  seen  that  the  gas  which  is  emitted  from  the  flame 
at  low  pressure  is  compressed  by  subsequent  pressure  waves  and 
raised  to  temperatures  higher  than  the  flame  temperature.  At  the 
same  time,  gas  emitted  at  low  pressure  is  later  cooled  below  the 
fleune  temperature  by  rarefactions.  Accordingly,  at  any  moment. 


McClure  assumes  that  a  particle  of  gas  traversing  the 
flame  is  subject  to  an  adiabatic  compression,  so  thati 

sLp  . 

eV*' 

which  may  be  substituted  into  equation  3  to  glve.i 


Then,  writing  Yp/v  «  l/f  —  dm/ifi/dp/p  and  ■  dp/p  ,  leads 
to  equation  6. 
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there  ^ould  exiet  in  the  duct  wevee  of  tenperature  which  ahould 
be  observable  by  their  luminosity.  An  idealised  representation 
of  the  appearance  of  the  gases  in  the  duct  is  shown  in  figure  4. 

An  alternate  description  of  the  condition  of  the  gases 
in  the  duct  may  be  made  through  the  introduction  of  the  entropy 


%nive  concept. 

for  a  polytropic  gast 

8  «  8^  4-  C  In  T/f®  -  R  In  P/f®  (7) 

P 

where  the  syndt>ols  have  their  usual  meanings.  If  it  is  assumed 
that  the  particles  are  all  emitted  at  the  same  flame  temperatw 
then  at  emission t 

8^  >  Constant  f  -  R  In  (6) 

but  in  an  oscillating  pressure  field: 

Pe  -  Pavg  (  1  +  f  ainoJ  t)  (9) 

substituting:  •  Constant  •  R  In  (1  £  sin  cj  t)  (10) 

Since  In  (1  ^  x)  >  x,  for  x  ^  1 

Se  ■  Constant  -  R  ^  sino^  t 

•  R  £■  sin  (  ^  t  *  TT  )  +  Constant 


Thus,  by  this  simple  relationshifik  the  pressure  and 
entropy  are  expected  to  be  180^  out  of  phase.  In  other  words, 
the  high  entropy  waves  are  emitted  at  pressure  minima.  The 
foregoing  derivation  assiuned  that  time  of  combustion  of  a  given 
particle  was  short  compared  to  the  time  of  oscillation.  If  the 
frequency  of  oscillation  is  very  high,  this  simple  derivation  doe's 
not  apply.  The  more  complete  derivation  is  being  prepared  for 
publication  as  a  journal  article. 

Experiments  are  now  being  planned  to  test  for  the  occur¬ 
rence  of  entropy  waves  using  the  previously  described  oscillator  (2). 
The  use  of  back-lighting  should  help  to  resolve  the  possibility 
of  "smoke",  since  the  transmission  of  light  through  the  combustion 
products  would  vary  with  "smoke"  density. 


An  even  more  significant  experiment  may  be  employed  to 
.test  the  entropy  wave  hypothesis.  A  schematic  of  the  equipment 
is  shown  in  Figure  5  .  The  experiment  is  effectively  a  time-resolved 
sodium  "D"  line  reversal  method  for  determining  gas  temperature, 
and  is  similar  to  a  method  used  to  study  the  temperature  history 
behind  a  shock  wave  (3)!.  For  this  purpose,  small  amounts  of  a 
sodium  salt  could  be  incorporated  in  the  propellant.  The  source 
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temperature  would  be  varied  to  obtain  an  indication  of  the  gas 
tcroporaturo.  The  dioico  of  a  cource  will  bo  dictated  by  the  flame 
temperature  of  the  propellant  under  study.  A  consideration  of 
the  entropy  wave  concept  and  adiabatic  coapreosion  indicates  that 
a  20%  temperature  variation  is  to  bo  e^ipected.  If  the  initial 
results  prove  presaising,  oUior  eninsrimonts  might  be  conducted. 
Among  these  are  first,  a  variation  of  the  source  temperature 
during  a  run  in  order  to  obtain  a  better  measurement  of  the  abso¬ 
lute  temperature  maonitudo.  A  second  teclvniquo  involves  the  use 
of  interference  filters  throbgh  whidt  the  gas  luminosity  could 
bo  photographed  so  that  the  temperature  variations  %rould  be 
accentuated  by  the  intensity  of  t)to  sodium  emission  recorded  on 
the  film. 

The  pressure  relationship  can  bo  obtained  directly  by 
the  use  of  a  gauge  which  can  be  mounted  in  the  tost  section  in 
the  same  plane  as  the  three  view  windows. 


V.  PARTICLE  TRACK  fIBTHOD 

Observations  of  velocity  profiles  by  means  of  luminous 
or  illuminated  particles  have  been  used  by  some  experimentalists 
to  trace  temperature  profiles  in  laminar  flames  (4).  The  shapes 
of  the  flow  pattern  around  various  wings  in  low-speed  wind  tunnels 
have  been  studied  with  luminous  tracers  ejected  into  the  gas 
stream (5).  Tliis  type  of  technique  has  been  considered  as  a  pos¬ 
sible  method  of  determining  the  response  function  of  the  burning 
propellant  surface. 

The  use  of  the  track  method  consists  of  determining  the 
velocity  of  the  tracer  as  a  function  of  time  and  position  by  eitheri 
(a)  following  individual  particles  through  a  number  of  frames  from 
a  movie  film  and  determining  velocity  from  the  displacement-time 
history,  (b)  illuminating  the  flov/  field  stroboscopically  and 
determining  the  different  velocities  from  the  length  of  the  lines 
\diich  result  on  the  film,  or  (c)  follov^ing  the  progress  of  the 
tracer  by  a  continuous  writing  streak  or  strip  camera,  and  deter¬ 
mining  velocities  from  measurements  of  the  slope  of  the  curve. 

A  knowledge  of  the  velocity  and  pressure  as  a  function  of  time 
gives  either  the  admittance  directly  or,  together  with  the  continuity 
of  mass  equation,  gives  a  knowledge  of  the  mass  flow  rate  as  a 
function  of  time.  Hence,  the  reaction  of  the  propellant  to  pres¬ 
sure  fluctuations  may  be  readily  determined. 

. .  The  feasibility  of  using  the  particle  track  method  will 

depend  upon  the  correctness  of  results  to  be  expected.  To 
investigate  the  fidelity  with  which  the  particles  follow  the  gas 
velocity,  consider  the  classic  Stokes  equation  for  the  equation 
of  motion  of  a  spherical  particle  in  a  gas  flowing  with  an  imposed 
sinusoidal  velocity  variation. 

:  0(  (y  ^  V  Oz) 


9. 


where  V  »  average  gao  velocity 

V  m  amplitude  of  gas  velocity  variation 
or  * 

yu  m  gas  viscosity 
D  «  particle  diameter 
«  particle  density 

The  solution  of  the  differential  equation  IIZ)  may 
bo  witteni  a  t/ 

/Jta  *  J  7==^  ^  (13) 

of*-  *u}^ 


where I 


('^1 


The  last  term  in  equation  (#3)  is  small  compared  with  the  first 
two  terms  since  the  magnitude  of  tt  is  such  that  the  ex¬ 
ponential  decays  very  rapidly  and  the  pre-exponential  terra  is 
less  than  imity.  . 

I  • 

Equation  (13)  indicates  that  three  factors  will  affect 
the  accuracy  with  which  observed  particle  tracks  will  reflect 
the  true  gas  velocity  history.  In  the  first  term  of  equation  («) ; 
the  expression  (  )  represents  the  extent  to  which  the 

particle  attains  the  steady  gas  velocity.  In  the  second  term, 
the  expression  (  1  +  to /v*-  represents  the  accuracy  with 
which  the  particle  track  will  indicate  the  amplitude  of  the  gas 
velocity  oscillations,  while  ^  represents  the  phase  lag 
between  the  gas  velocity  oscillation  and  the  particle  velocity 
oscillation.  The  effects  of  these  factors  may  be  considered 
separately. 


To  study  the  error  from  the  first  term,  the  time  and 
distance  required  for  a  particle  to  reach  9596  of  the  gas  velocity 
were  calculated.  Using  a  velocity  of  500  cra/sec.,  /O  .  2.8,  and 
«  .063  centipoise.  The  results  are  as  follows  i 


Particle  Diameter  (microns) 
5 

10 

20 


T'  (msecs.)  ^  (cm.) 
0.24  0.08 

0.95  0.32 

3.79  1.30 


It  should  be  mentioned,  at  this  point,  that  the  dis¬ 
tances  resulting  are  large  compared  to  the  combustion  zone  thick¬ 
ness  (.Olcra),  thereby  justifying  the  use  of  a  steady  flow,  rather 
than  a  flow  which  initially  undergoes  a  very  steep  acceleration  to 
500  cm/sec.  in  .01cm.  This  assumption  greatly  simplifies  the 
equations  and  resulting  calculation. 
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To  study  tho  error  from  tho  second  term,  calculations 
were  made  to  determine  tho  magnitude  of  the  phase  lag  as  well  as 
the  amplitude  of  rft  ccnqpared  to  V  •  The  results  aro  as  fellowst 


0  (microns) 

Phase 

Lag 

Relativo  Amplitude 

100  qps 

SOOeps 

100  cps 

SOOqps 

5 

2.80 

13.90 

1.00 

.98 

10 

11.20 

44.7® 

.98 

.71 

20 

38.40 

75.90 

.78 

.45 

Zf  m  content  oursolvos  with  a  90X  relative  anqplitude, 
then  the  criterion  that  0.5*  can  be  imposed.  This  criterion 
leads  to  a  maximum  limit  in  phase  error,  and  also  leads  to  the 
particle,  attaining  95%  of  steady  gas  flow  velocity  in  less  than*% 
period.  Substituting,  tho  criterion  loads  to  the  requirement  that! 

(D  in  microns) 

The  results  are  plotted  in  Figure  7 .  The  effect  of  tho  density 
of  the  particle  should  be  noted.  As  a  side  note,  the  imposition 
of  95%  relative  velocity  amplitude  leads  to  w/er  <  e.f ,  giving  a 
phase  error  of  about  and  catching  up  times  on  the  oi^er  of 
1/20  period.  This  criterion  drops  the  maximum  frequency  for  a 
given  particle  by  a  factor  of  5. 

Consideration  of  the  calculations,  reveals  rapidly  that 
very  small  particles  must  be  made  visible  for  the  method  to  have 
any  feasibility.  As  a  result,  good  lighting  and  excellent  optics 
will  be  required.  First  attempts  will  be  made  with  the  set-up 
shown  in  Figure  6  .  A  flash  bulb  might  be  used  because  25msecs. 
would  give  10  periods  at  500  cps*  and  allow  5  separate  "families* 
of  particles  to  be  followed  from  the  surface  through  the  viewing 
section.  The  use  of  a  chopper  to  give  periodic  illumination  malMS 
velocity  determinations  far  simpler,  but  steady  light  could  be 
used  to  give  a  complete  history  of  a  particle,  even  though 
absolute  measurements  are  harder.  The  focusing  can  be  varied 
depending  on  the  light  requirements.  In  the  drawing,  a  "pencil* 
of  light  is  used,  so  that  fewer  particles  will  be  under  observation, 
leading  to  greater  clarity  in  interpretation.  Focusing  on  a 
point  as  opposed  to  a  line  might  give  better  results. 

Experiments  will  initially  be  made  using  well-screened 
MgO,  since  this  oxide  has  a  melting  point  (3100^  K)  in  excess  of 
flame  temperatures  and  should  thus  keep  its  shape  well  in  the 
combustion  gases. 

VI.  CONCLUDING  REMARKS i 

The  general  objectives  of  the  experimental  program 
during  the  next  quarter  are  as  follows i 
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1.  The  operating  parameters  of  the  oscillator  driver  will  be 
investigated. 

2.  The  combustion  of  propellants  of  various  types  under  conditiens 
of  oscillating  pressure  will  be  studied  by  direct  photographs. 

3.  Experiments  to  measure  gas  temperature  and  gas  particle 
velocities  will  be  instituted. 


# 
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TABLE  I 


BxperlMBtal  conditions  under  which  film  strip  was  ss^osodi 


Moan  ^anbsr  prsssurs 
Frsssurs  oscillation  anplituds 
Msan  burning  rats 
Vrsqusncy  of  oscillation 
Csmsra  speed 
Propellant 

Theoretical  flame  temperature 


-  500  psi 
±  125  psi 
-0.5  in./sec. 

-  500  qps. 

->  5000  frames/sec. 

-  composite  double-base, 
containing  2%  Mg 
(through  400  mesh) 

-  3100®  K 


i 


1 


In  this  figure  the  pressure  and  temperature  of  gas  packets  emitted 
from  the  flame  are  plotted  as  a  function  of  time.  Particles  emit¬ 
ted  at  various  times  during  the  pressure  cycle  all  have  the  same 
temperature.  However,  when  they  enter  the  duct  their  different 
entropies  are  manifested  as  temperature  differences.  Thus,  gas 
particle  1,  emitted  at  high  pressure  and  low  entropy,  senses  an 
expansion  and  cooling  to  some  lower  temperature,  and  its  temperature 
oscillates  between  this  temperature  and  the  flame  temperature. 
Particle  3,  emitted  at  low  pressure  and  thus  at  high  entropy,  is 
subsequently  compressed  and  heated  to  some  higher  temperature.  The 
temperature  of  this  particle  emitted  at  low  pressure  oscillates 
between  this  higher  temperature  and  the  flame  temperature. 
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FIGURE  7 


In  this  figure  are  plotted  the  maximura  sizes  of 
particles  which  will  faithfully  (within  90%)  of 
true  amplitude)  follow  a  velocity  oscillation  of 
a  given  frequency. 
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